WUB/09-05 



Mean-Field Gauge Interactions in Five Dimensions I. 

The Torus 



Nikos Irges and Francesco Knechtli 

Department of Physics, Bergische Universitdt Wuppertal 
Gaussstr. 20, D-42119 Wuppertal, Germany 
e-mail: irges, knechtli@physik.uni-wuppertal.de 



Abstract 

We consider the lattice regularization of a five dimensional SU{2) gauge theory with 
periodic boundary conditions. We determine a consistent mean-field background and 
perform computations of various observables originating from fluctuations around this 
background. Our aim is to extract the properties of the system in regimes of its phase 
diagram where it seems to be in a dimensionally reduced state. Within the mean-field 
theory we establish the existence of a second order phase transition at finite value of 
the gauge coupling for anisotropy parameter less than one, where there is evidence for 
dimensional reduction. 
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1 Introduction 



It is possible that our world has more than four space-time dimensions. There are different 
ways that extra dimensions could leave a trace depending, among other things, on which 
of the fundamental forces we choose to look at. We will concentrate here on the gauge 
interactions. Leaving to the side for the moment the fermionic sector we present here 
our investigation of the dynamics of extra dimensional pure gauge theories with focus on 
dimensional reduction. More specifically, we would like to propose a scheme where this 
phenomenon can be analyzed analytically, far from the perturbative regime. 

The phase diagram of a five dimensional SU (N) gauge theory in infinite four dimen- 
sional spatial volume is parametrized by the two dimensionless couplings (3 = 2N/ (g^A) 
and N5 = 27iRA, with the five dimensional bare coupling, R the physical length 
parametrizing the size of the fifth dimension and A a cut-off q The domain of weak 
coupling perturbation theory is the vicinity of the "trivial" point P, ^ 00 of the phase 
diagram which is approached as one removes the cutoff. From the formal point of view, 
up to now most analytical investigations of higher dimensional theories have been carried 
out in this domain, which however comes with two caveats: as one approaches the fixed 
point, physics is governed by the triviality of the coupling and as one tries to enter in 
the interior of the phase diagram cut-off effects dominate. Far from the trivial point, in 
the small /3 regime, the system finds itself in a confined phase below a critical value (3c of 
order one. The regime Pc < P << 00 where neither perturbation theory nor the strong 
coupling expansion is useful, up to now, is analytically basically unexplored. 

The success of the Standard Model (SM) in explaining experimental data requires, 
after introducing a new ingredient, a natural way to sufficiently hide it. As far as higher 
dimensional theories are concerned this would mean either that the extra dimensions are 
very small or that there is a four dimensional slice in the higher dimensional space where at 
least part of the physics is localized. From the phenomenological point of view the former 
situation has been also investigated by standard methods, mainly with a combination 
of Kaluza-Klein theory and weak coupling perturbation theory In this approach the 
size of the extra dimension(s) can be tuned to small values and the delicate issue is 
to make it small enough so that it does not lead to a contradiction with the data but 
large enough so that it can have a sizable effect. In the localization scheme, the extra 
dimension can be, in principle, as large as any other dimension (or even larger) but a 
dynamical mechanism is necessary to implement it. Unlike for gravity, for gauge fields 
a classical localization mechanism is not known and perturbation theory does not seem 
to help in this respect. Consequently in most model building approaches when a certain 

^ Later we will be interested in anisotropic spaces where a third parameter, 7, will appear. 
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field needs to be localized (for some phenomenological reason) one just assumes that it 
is. Nevertheless, a few non-perturbative quantum mechanisms of gauge field localization 
were invented in the past [2l[3]. In these, lattice Monte Carlo simulations have been proven 
a crucial tool. Despite all efforts, from the five dimensional point of view a quantum field 
theoretical approach to the localization of non-abelian gauge fields has been elusive. 

This state of affairs calls for a complementary tool to the perturbative computations 
and the numerical simulations, one which would enable us to probe the system analytically 
up to its phase transition. The method known to work in the domain of /3 of order one 
is the mean-field approximation |11|5]. This approach has been used before to locate the 
critical value of f3c for various gauge theories on toroidal geometries (see |1] and references 
therein) and in [H] it was used to locate /3c for an orbifold geometry. In the 80's there was 
a considerable effort to define an expansion in the fluctuations around the saddle-point 
(or zero'th order) approximation [7]. The free energy was indeed computed to a high 
order in order to locate any phase transitions very precisely and hopefully predict their 
order. However, extensive observables other than the free energy were never computed 
sytematically especially for theories of large dimensionality. 

In fact, it is believed that the expansion around the mean- field background becomes 
a better and better approximation to the non-perturbative behavior of the system as 
the number of space-time dimensions increases so it seems that it is just tailored for our 
purposes. In addition, we need a regulator that maintains a finite cutoff while preserving 
gauge invariance. The computational scheme will be therefore a lattice regularization 
with lattice spacing a = 1/A in the mean-field approximation to zero'th order. At higher 
orders we will perturb away from the mean-field background by corrections that indeed go 
as one over the number of space-time dimensions. For the basic definitions and conditions 
for a consistent formulation of a lattice non-abelian gauge theory in five dimensions as 
well as some related Monte Carlo results see [8],[9]. 

In section 2 we describe the mean-field method applied to gauge theories and show 
how to compute local and global observables using fluctuations around the background. 
In section 3 we turn to the definition of lattice gauge theories around the mean-field 
background and concentrate on SU{2). In section 4 we derive the free energy, the static 
potential and the mass formulae for the SU{2) lattice gauge theory in five dimensions. 
In section 5 we discuss numerical results for these quantities and determine the phase 
diagram. In section 6 we give our conclusions. 
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2 Gauge theories in and around the mean-field background 



2.1 The zero'th order approximation 

The partition function of an SU (N) gauge theory on a Euclidean lattice is 

Dt/e-^«[^l , (2.1) 



where is the Wilson plaquette action defined in terms of oriented plaquettes U{p). 

A plaquette is the product of four links around an elementary square on the lattice. Links 
are generically denoted by U. The mean-field computation starts by inserting an identity 
in the path integral, using the Fourier representation of the 6 function 

6{f{x)) = ^|l^e-(^)^(^). (2.2) 

This allows to integrate out the SU{N) variables U in favor of a set of unconstrained 
variables V and the set of Lagrange multipliers H that implement the constraint U = V: 

Z = j J^iU j T)V 5{V - U)e-^^^^^ (2.3) 

= JdU JdV J Di7e(^/^)^^*'-{^(^-^»e-^«[^]. (2.4) 

An effective mean-field action Sc{{[V,H] is defined using 

e-iH) = j j^u^(l/N)Retr{UH} ^2.5) 

in terms of which we can express conveniently the partition function as 

Z = J^V J Di/e-^^ffl^'^l , ^eff = Sg[V] + u{H) + (l/A^)Retr{W} . (2.6) 

The mean-field or zero'th order approximation amounts to finding the minimum of the 
effective action when 

H — >H1, V -^Vl, Ses[V, H] =minimal . (2.7) 

Taking derivatives of eq. (12. 6p with respect to V and H one obtains the equations 

V = 



OH 

H = - 



H 

dV 



(2.8) 

V 
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which determine the mean-field background. It is a set of non-hnear coupled algebraic 
equations that can be solved numerically by an iterative method. 

The free energy per lattice site is defined as 



F = --^ln(Z), 

where M is the number of lattice sites. At 0th order it is simply 



(2.9) 



(2.10) 



2.2 First order corrections 

Gaussian fiuctuations are defined by setting 

H = H + h and V = V + v 
and Taylor expanding the effective action to second order in the derivatives § 



Scs = See[V, H] + 



1 / 5^S,s 



6H^ 



6H6V 



hv + 



cff 



V,H 



We define the quadratic pieces of the propagator 



eflf 



6H^ 

'eff 



V,H 



SVSH 

'eflf 



6V^ 



2 T^i'"'") T T^(vv) 

V = ViKl- 'vj = V 't 



V,H 



(2.11) 



(2.12) 



(2.13) 
(2.14) 
(2.15) 



in terms of which we can express the part of the effective action quadratic in the fluctua- 
tions as 

5(2) [^;, /i] = 1 [h^K^f^'^^h + 2v^K^''''^h + v^K^^^^v) . (2.16) 
The integral of the fluctuations 



Dv / D/ie-^'''[^'^] 



is a Gaussian integral and it can be easily performed to give 

_ (27r)|fc|/2(27r)l^'l/2 



(2.17) 



(2.1^ 



^The first order correction to the mean- field approximation in the fluctuations is second order with 
respect to the derivative expansion. 
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with \h\ and \v\ the dimensionahties of the h and v parameters respectively and 

K = + _ (2.19) 

Therefore to first order 



The typical quantity of interest is 



Df/C[f/]e-^G[t/] 



(2.20) 

(2.21) 
(2.22) 



with O a gauge invariant operator. Note that the 5-function constraint in Eq. (12. 3p 
has been used to replace 0[U] with 0[V]. The purpose of this step is to replace a 
constrained SU{N) path integral by a path integral over unconstrained matrix valued 
complex numbers. The number of degrees of freedom apparently increases but in fact the 
relevant constraint is encoded in K^'^'^\ 



The observable has a Taylor expansion of the form 



0[V] = 0[V] + ^ 



1 s^o 



V 



The combined expansion of observable and action to the same order is 



(O) 



0[V] 



2SV^ 











z J 




V 





V 

V 



The link two point function can be integrated to 



< ViVj >= - I Dv I DhviVie 



-5(2) [v,h] 



-1^ 



l%] 1 



(2.23) 

(2.24) 
(2.25) 

(2.26) 



where K is the matrix defined in Eq. f l2.19p . The indices i and j denote collectively the 
indices on which the links depend. The physical observable expectation value to first 
order can be then expressed as 



<0>= 0[V] + ^tr 



At first order we define the free energy from Eq. (12.201) as 









V ) 



= - ln(z) = + 4? In 



(2.27) 



(2.28) 



where we have dropped an irrelevant additive constant and App is the Faddeev-Popov 
determinant that appears after fixing the gauge. 
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2.2.1 Phase transitions 



At zero'th order in the confined phase, the free energy F has a minimum at if = (that 
typically implies = 0) which dominates over any other minima located at any other 
finite values of H. As (3 increases, a second local minimum develops and the critical value 
of P, which indicates when the system enters in the deconfined phase, is reached when 
this second local minimum becomes of equal height as the minimum at if = 0. In the 
deconfined phase, the local minimum at ii 7^ turns into a global minimum. Typically, 
between the two minima the free energy develops a local maximum. 

An alternative way to locate the phase transition is to look at the zero'th order mean- 
field conditions Eqs. (12. 8p and define the deconfined phase wherever the numerical it- 
erative method converges to a non-trivial solution. In general however we would like to 
emphasize that this is a safe method only as long as one is sure that the solution found 
corresponds to a global minimum of the free energy. Also we point out that the numerical 
value of f3c depends quite a bit on wether we fix the gauge or not (this is an option at 
zero'th order), wether we use the iterative method or the direct minimization of the free 
energy, but not its existence. The numerical value has no physical meaning as far as we 
can tell. The only thing that has a physical meaning is if we are "near" or "far" from the 
phase transition and whether the phase we are sitting in is the global minum of the free 
energy or not. 

To first order it is however necessary to fix the gauge. In gauges which require ghosts, 
a divergence in the free energy develops due to the presence of the Faddeev-Popov deter- 
minant that vanishes when the background is taken to zero. If the divergence is regulated 
(i.e. dropped) then the phase transition can be still located at the value of f3 below which 
the global minimum is at = 0. If it is kept then the global minimum of the free energy 
has a discontinuous jump from a larger value to a much smaller but non-zero value. The 
value of (3 at which this jump takes place signals the phase transition. 

Our last comment on the issue of the phase transition concerns its degree. In [10] it 
was established and in |3] confirmed that on the isotropic lattice it is a strong, first order 
phase transition. 

2.2.2 Extracting masses at first order 

Let us take a generic, gauge invariant, time dependent operator or "observable" 0{t) and 
its 2-point function O{to + t)0{to). In order to extract the mass associated with 0{t), we 
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need the connected correlator 



C{t) = < O{to + t)C(to) > - < O{to + t) >< O{to) > 

= C(o)(t) + C(i)(t) + ■■■ (2.29) 



where the C^^\t) and C'^^^t) can be identified from 



<0{to + t)0{to)> = 0W(to + ^)0^°H^o) + ^tr{ ^'^^^^°+^^^^^^°^^ X-^} 



H 

(2.30) 



Since C'^°^(t) = and because time independent contributions drop out from C^^\t), at 
first order we have 

2 [ d-^v J 2 I o^'i' I 

where the notation 5*^^'^^ means one derivative acting on each of the O{to + 1) and O(to) 
operators. In the second part of the equation tilded quantities are the Fourier transforms 
of the corresponding untilded quantities. 

To extract the scalar mass we observe that a gauge invariant correlator admits an 
expansion in terms of the energy eigenvalues of the states it contains as 

C(t) = J]c,e-^^*, (2.32) 

where Eq = m , Ei = m* , ■ ■ ■ and therefore we have that 

m ~ lim In ^ \ . (2.33) 

t^oo C(l)(t-1) ^ ^ 



2.2.3 The static potential 

The static potential can be obtained from a Wilson loop Ow extending in the time 
direction and in one spatial dimension. In the infinite time limit the expectation value of 
the loop is related to the static potential as 

t^oo: e"^* ^<Ow> ■ (2.34) 

This immediately implies that in the mean-field background the potential is constant 

Ow[V] = Nvl^'^+'\ (2.35) 

To obtain the first order correction, one must compute a gauge boson exchange between 
the two time-hke legs of the loop plus the self energy and tadpole diagrams. All other 
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possible exchanges of gauge bosons vanish in the infinite time hmit. The static potential 
is the quantity that will decide whether the system is in a dimensionally reduced state or 
not. 

The static potential can be extracted to first order form 

V=-\im^- log iO^iV]) - hm (2.36) 

t^oot ^ ^ t^oct Ow[V] 

with the correction given by the sum of the exchange, self energy and tadpole diagrams. 
2.3 Second order corrections 

We will encounter cases where a physical observable (like a gauge boson mass for example) 
due to the choice of the operator that represents it, is identically vanishing to first order. 
We have to know therefore how to generalize our computations to second order in the 
expansion around the mean-field background. 

To second order we have the expansion 

S^, = S^,[V, H] + - ^-^h + 2^hv + J 

(2.37) 

The cross terms in the cubic and quartic terms vanish because of the special form of S'eff. 
We have a similar expansion for the observable 

^^^^^^ SO 16^0 2 IS'^O , 1 5^0 4 , , 

OIV] = OlV]^-v+-^v^ + -^v^ + -^v' + - (2.38) 

Straightforward algebra leads to the result that the tadpole-free contributions to the 
observable are 

< O >= Om + \ 1 DhDv {^^^v^ + e-^'"|.,ft| (2.39) 

where the second order correction involves the integral 

^ j Dv j Dh {viVjViVm)e~'^''^\v,h] = 

{K-%{K-')i^ + {K-%{K-'),^ + {K-\m{K-')ji (2.40) 

The contribution leads to a time independent tadpole contribution to the correla- 

tor Eq. (12.291) and can be dropped. We can finally express the physical expectation value 
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as 



(2.41) 



To extract the masses to second order is straightforward. Again, all time independent self 
energies cancel from connected correlators and at the end the mass is obtained from 

" - ,'™ C<')(« - 1) + CP)(« - 1)' ^ ' 

with 



i,j,l,m ^ ^ ^i^'" 



(2.43) 



2.4 Gauge dependence and zero modes 

Gauge dependence of the results turns out to be the most controversial issue in computing 
fluctuations around the mean-field background. It can be shown [H] that for any gauge 
invariant polynomial P{U) of the links it is true that 

^ j DUe-^'''P{U) = ^ J ^U5{Ui,g)e-^'^P{U) (2.44) 

with Sw the gauge invariant (for example Wilson plaquette) action, and g an SU{N) 
group element. The inserted term 6{Ui,g) allows the link Ui to be set equal to g. This 
implies that one can set any particular link to a chosen element which is typically taken to 
be the unit element. Iterating this construction one can set to unity any number of links, 
as long as there are no closed loops forming in the process. The argument then is that 
since the lattice action, the measure and the operator "observable" P{U) are explicitly 
gauge invariant, then any correlator is also gauge invariant. This allows one to compute 
the correlator in a fixed gauge and the previous argument guarantees that this will not 
affect the gauge invariant result, as long as one can compute exactly the left hand side of 
eq. ( ]2.44p analytically. This is however not possible in general and one either computes 
the correlator via a Monte Carlo simulation or analytically in some approximation. 

Higher dimensional gauge theories offer an excellent opportunity to combine success- 
fully the Monte Carlo and mean-field methods because of the expectation that the expan- 
sion around the mean-field background becomes a better approximation as the number 
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of space time dimensions increases Thus, in a five dimensional gauge tlieory one can 
compute the same quantities in two different ways, one numerical and gauge invariant 
and one analytical but gauge dependent and if these two agree to a satisfactory degree 
(as one would expect from general arguments) one can be confident that the calculations 
are not out of control. By "satisfactory degree" we mean that both methods should de- 
scribe the same physics qualitatively (this is a necessary condition) and they must be in 
a good quantitative agreement. By good quantitative agreement we mean that of course 
the numbers should not be off by order of magnitudes and if the two methods agree on 
the physics, some quantitative disagreement can be perhaps tolerated. The next issue is 
which gauge to choose. 

Taking maximal advantage of the gauge fixing possibility amounts to fixing the links 
to unity on a spanning (or maximal) tree. The disadvantage of this "maximal" gauge is 
that it is hard to perform explicit calculations with it for non-Abelian groups. Another 
simple choice is that of the axial gauge where one fixes to unity links along only a specific 
direction. An advantage of this gauge is that it is by construction ghost free. A great dis- 
advantage is that it leaves part of the local gauge invariance unfixed and as a consequence 
most observables are plagued by zero modes of the propagator which correspond to local 
gauge transformations that are independent of the gauge fixed coordinate. In momentum 
space they appear whenever the momentum pj = 0, where j is the gauge fixed direction. 
If one restricts the analysis to observables that do not see these zero modes (i.e. if they 
couple to sectors of the propagator with non-zero eigenvalues only), the axial gauge is 
actually fine. If the local zero modes do affect the observables, one must fix the residual 
gauge and then ghosts appear. A further subtle issue of the axial gauge concerns the 
j-coordinate dependence of the various quantities in the axial gauge on a finite lattice. In 
order that the axial gauge is performed, one has to assume that the gauge fixed direction 
is infinite, however one can compute in practice the propagator and observables only at 
a finite extent of the j direction. In order to minimize the ambiguities originating from 
this, one should take this direction as large as possible. 

Finally, there are more conventional gauges such as the Lorentz gauge which is the one 
we will use here. It involves ghosts and one must take account of them. All gauge invariant 
correlators evaluated in the zero'th order background turn out to be gauge independent. 
In some cases, like for the static potential and the scalar field's mass, we will be able to 
show this analytically and in others, like for the vector field's mass, we will be able to 
show it only numerically. The only quantity which will have some gauge dependence is 
the first order correction to the free energy. We do not fix the gauge at zero'th order and 
we will make sure that the gauge dependence of the first order correction results only in 
quantitative changes (i.e. it moves around a bit the minimum) but does not change the 
qualitative properties of the system, for moderate varying of the gauge fixing parameter. 
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This of course means that numbers originating from the minimization of the free energy 
should not be given any physical meaning. If this had to be the case, like in a next to 
leading order calculation, a more careful treatement would be necessary. Since we will 
not use the first order free energy for anything else other than showing stability, we will 
not need to carry out any special regularization. 

As already stated in [4J, in any gauge with periodic boundary conditions there are 
additional zero modes, the global zero modes (torons) which appear when all lattice mo- 
menta vanish simultaneously. These zero modes are present and they must be taken into 
account as well [12]. In our case, they always show up as a zero of the inverse propagator 
at the same time with a zero in either the observable or in the Fadeev-Popov determinant. 
The result will be always a zero over zero contribution which can be regularized and a fi- 
nite part can be extracted. This finite piece however is volume suppressed and contributes 
a negligible amount. Our regularization of the toron contributions will be thus to simply 
drop them. The only other regularization we do is drop an overall constant contribution 
from the free energy and which has no physical meaning. 



3 The lattice model and its observables 
3.1 The action 

The discretized version of the torus is defined on a five-dimensional Euclidean lattice 
with lattice spacing a. The points are labeled by integer coordinates n = {n^/}. The 
dimensionless lengths of the lattice are L = l/a in the spatial directions (M = /i = 1, 2, 3), 
T in the time-like direction (M = /i = 0) and = 271 R/ a in the extra dimension (M = 
5). Periodic boundary conditions are used in all directions. The dimensionless length of 
a direction on the lattice is just the number of points in that direction. Occasionally, we 
will use the notation A/" for the total number of points in the lattice. The convention 
for the fifth coordinate labels is that = 0, ■ ■ ■ , — 1 . The lattice momenta^] are 
corresp ondingly 

PM = -j—k, /c = 0, 1, . . . ,Lm - 1 (3.1) 

and Lm is the dimensionless length of the lattice in direction M. The gauge field is the 
set of link variables {U{n, M) G SU{N)} and its action is taken to be the Wilson action 

Sw[U] = ^J2^'i^-UiP)}^ (3.2) 
p 

^ Unless otherwise stated, dimensionless coordinates and momenta are used. 
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where the sum runs over all oriented plaquettes U{p). A plaquette at point n in directions 
M and K is defined as the product 

U{p) = ti ^^U{n,M)U{n + M,K)U\n + k,M)U^{n,K)^ (3.3) 

of links. A gauge transformation Q acts on a link as 

f/(n, M) — > n{n)U{n, M)n^{n + M) (3.4) 

The Wilson action reproduces the correct naive continuum gauge action of the fields 
Afj, B = 1,...,N'^ — 1 of the five- dimensional gauge potential Am defined through 
U{n,M) = exp{aAM{n)}. 

We will now apply the general formalism we described to a specific example: an SU (2) 
lattice gauge theory in 5 dimensions with fully periodic boundary conditions. 

3.2 Ghosts 

A general 2x2 complex matrix can be represented by 

3 

Vo + i'^VAa'^, Vo,VAeC. (3.5) 

A=l 

We choose the Lorentz gauge, which amounts to adding the gauge fixing term 

'5gf = ^EE[/^H]' (3-6) 

A=l n 

to the lattice action. We will now derive the corresponding Faddeev-Popov term. Notice 
that we are fixing the link fiuctuations, not just the gauge field fiuctuations. Also, we do 
not fix the link fiuctuations along the Hermitian component vq. Consistently, it will turn 
out that the propagator for vq does not have poles. A similar statement is implied in |1] 
for SU{N) in general. 

In order to compute the Faddeev-Popov determinant we have to compute the variation 
of the gauge fixing term 

fA{n) = ^K(n, M) - VA{n - M, M)] (3.7) 

M 

under infinitesimal gauge transformations. We parametrize gauge transformation func- 
tions as n(n) = e*'^'*(")°"^ = l-\-iuj^{n)cr^ where the second equality holds for infinitesimal 
transformations. Then, the gauge transformation rules of the components of a link are 

6vo{n,M) = v^{n,M) (u^{n + M) -io'^in)^ 
5vc{n,M) = -vo{n, M) (^u^ {n + M) - uj^\n)^ 

- e^^S^(n,M) (^a;^(n + M)+a;^(n)) (3.8) 



13 



and the transformation they induce on the gauge fixing function is 

Sfcin) = [Svc{n,M)-6vc{n-M,M)j = J^^Cn.AmU^^M (3-9) 



in,A 



with the ghost kernel defined as M.cn;Am = ■^Cn'-Am 



■^^'■Cn;Am 



S^^ ( voin, M) + fo(n - M, M) 



— e 



ABC 



+ 6. 



m,n+M 



m,n—M 



VB{n, M) -VB{n- M, M) 



[-5^^t;o(n, M) - e^^S^(n, M)] 



5^\Q{n -M,M) + e^^^v^{n - M, M) 



ABC„.Bi 



The ghost action is then 



^FP = Yc^{n)MAn;BmC^i 



m) 



(3.10) 



(3.11) 



The quadratic part of this defines the Faddeev-Popov kernel which in a general mean-field 
background V = V{n, M) ■ 1 is 

M. An:,Brn{y) = 

^^'^ \<rn (F(n, M) + Vin - M, M)) - 5^ F(n, M) - 5^ F(n - M, M) 

(3.12) 



Note that the gauge boson-gauge boson-ghost vertex is proportional to V. In a con- 
stant background Eq. fl3.12p is the lattice version of the Laplace operator, while in a 
non-constant background its appropriate generalization. Note also that 6vo{n)\y = is 
consistent with not gauge fixing vq. 



3.3 Observables 

We first define the auxiliary "lines'" 

715 — 1 



/(n5)( 



and 



ms=0 



to+t-1 



5) = J]^ ^VqI + Va {{mo, m,m5); 5) 

1715=0 



(3.13) 



to+t-l 



l^^\to,m,m5) = Yl U{{mo,m,m5);0) = JJ IvqI + iimo,m,m5);0) a"') , (3.14) 



mo=to 



mo=to 
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Figure 1: Contributions to the static potential on the torus: from left to right, gauge 
boson exchange, self energy and tadpole. 

where in the second parts of the equations we inserted the mean-field parametrization 
and we have introduced the matrices 

a° = {l, ia^}, a° = {l, -ia^}, a = 0,A (3.15) 

The line along the extra dimension can start from any spatial point m but will be always 
taken to start from = while the temporal line can start from any to,m,m5. The 
superscript indicates the extent of the line. Finally, we will be using the label mo or t,t', t" 
etc. for the temporal components according to convenience. 

The static potential will be computed using the averaged version over the 4d starting 
position (to, rh) 

C?r"^^ = ^E^!^"n^o,m) (3.16) 

to,m 

of the Wilson loop observable 

0%-'''\to,m) = ii[l^'\to,m,Q)l^'^^\t + to,m)l^^^ (3.17) 

Similarly to /^"^^ one can define spatial lines Z*^"*) along the /c = 1,2,3 directions as well. 
In fig. [1] the diagrams contributing to the static potential can be seen. 

For the masses we first define the Polyakov loop 

P(°)(t,m) = /(^s)(t,m) (3.18) 
in terms of which we define the operator 

$(°)(t,m) = P(°)(t,m) -P(°)^(t,m). (3.19) 
The scalar gauge invariant operator 

Co(t,m) = tr{p(°)(t,m)} (3.20) 
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is then used to determine the Scalar ("Higgs") □ observable, determined by the averaged 
over space and time location of the operator 

Onit) = ^J2j6j2 Oo{to,m')Oo{to + t,m") (3.21) 

to m' ,rn" 

Next, we construct the displaced Polyakov loop 

^'^(t, m) = U{{t, m, 0); k) $(°)^(t, m + k) U{{t, m, 0); k^ $(°)(t, m) (3.22) 
which is an operator with a vector and a gauge index. The gauge covariant quantity 

0,^(t,m)=tr{<)'^(t,m)} (3.23) 
will be used to define our Vector ( ' W ) mass observable, via the 2-point function 

W = ^ E E E E ^^(^0, rn') 6m Of{to + t, m") (3.24) 

to m' ,rn" A k,l 

where A = 1,2,3 is a sum over the gauge index and fc,/ = 1,2,3 a sum over the spa- 
tial Euclidean index. We have set the lattice spacing a = 1 for simplicity. It will be 
reintroduced when we discuss the results, where it has a physical meaning. 

The choice eq. fl3.20p gives a non-trivial result for the Higgs mass already at first order. 
A consequence of this choice is that since it is essentially a single gauge boson exchange 
diagram, see Fig. [21 there is no L-dependence built in it by construction. Therefore, by 
measuring the mass using this observable, we conveniently obtain its infinite L limit value. 
If desired, the second order correction could be added which will bring in an intrinsic L- 
dependence. Here, we will restrict ourselves to the first order expression. In Fig. [3] the 
second order contribution to the gauge boson mass is shown. 

These choices are of course not unique but not arbitrary either. Their general form is 
fixed by their transformation properties under charge conjugation (C), three dimensional 
coordinate inversion (P) and by their Lorentz (spin J) and gauge (isospin /) indices [13] . 
The basic quantity is the field which represents the extra dimensional component of 
the gauge field (which we have called Higgs for short). It is easy to see that its continuum 
limit is |9] 

^'^'^^ =4N5aA5 + 0{a^). (3.25) 

We will be interested in taking the continuum limit at finite R = ^r^. In this limit, 
is fixed. The continuum limit of is on the other hand [9] 

Ot = a tr{a^(Z}fc<l>(°)t)${o)| ^ (3.26) 

''We will call the scalar Higgs even though we will not be interested here in the phenomenon of 
spontaneous symmetry breaking. We have in mind the analogous observable on an orbifold geometry 
which in the presence of a vev for the scalar can trigger symmetry breaking |16] . 



16 



Figure 2: Contributions to the mass of the scalar at first order in the mean-field. The 
blobs are arbitrary lattice points, identified on the torus. 

to+t _ 




Figure 3: Contribution to the gauge boson mass at second order in the mean-field. Points 
at a fixed n^, at ^5 = and = are identified. 

where Dk^^*^"* = dk + [Ak, $''*^^] is the covariant derivative. 

The charge conjugation C leaves the coordinates invariant and acts on links as 

U* = a^Ua\ (3.27) 

while parity P acts on the coordinates as Pn = np = {no, — n, n^) and on the links as 

PU{n,k) = U\np-k,k), 
PU{n,0) = U{np,0), 

PU{n,5) = U{np,5). (3.28) 

By construction, Oq of Eq. fl3:20|) has C = P = 1 (scalar) and of Eq. fl3:23|) has 
C = 1,P = — 1 (vector). The isospin is defined by the transformation of the operators 
under global gauge transformations. For the definition of the spin of lattice operators we 
refer to [H]. It is easy to see that Oq has / = J = and Cj^ has I = J = 1. 
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4 Free energy, static potential and mass formulae in the SU(2) 
theory 



In this section we continue our discussion on the SU{2) theory. In expressions vahd 
though for general SU (N) we keep an exphcit N. 

4.1 At zero'th order 

The SU (2) hnk matrices have the representation 

3 

U{n,M) = uo{n,M)l + i^UA{n,M)a^ , Ua^R, UaUa = I , (4.1) 

A=l 

where a = 0, 1, 2, 3. The auxihary matrices V and H are arbitrary complex 2x2 matrices 
and can be represented by 4 complex numbers as 

V{n,M) = vo{n,M)l + i^VA{n,M)a'^ , Va e C , (4.2) 

A 

H{n,M) = ho{n,M)l-iJ2hA{n,M)a'^ , /i„ e C. (4.3) 

A 

A standard computation gives for the effective action u the result 



u 



(H(n, M)) = - In {^h{p)^ , p{n, M) = lY^{ReK{n,M)Y (4.4) 



with Ii the Bessel function of the / type and evidently it does not depend on the imaginary 
parts of the h^s,. The 5 function that we insert in the partition function is 



nr. A/f *^ 



n M a 

-|-|- -|-|- -|-|- /■+^°° dRe /i„dlm h 



-dRei;„dIm Vo 



2 

exp{-(Re/ia)[(Re^;a) - Ua] - {lmha){lmva)} ■ (4.5) 



Since u{H) does not depend on Im ha the integration over Imha will simply give 6{lmva) 
and eliminate these degrees of freedom from the computation so finally we can take 

Va^R and /i^ e M . (4.6) 



'S'eff = J] RetrV(n;M,i^) 



The effective action is therefore 

2 

+ E E I ^('"(^' ^)) + E ^-{n, M)va{n, M) ) , (4.7) 
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n M<K 



n M 



where V{n; M, K) is the plaquette product at position n in directions M and K. When 
we calculate corrections 



we set 

ho{n,M) — > ho + ho{n,M), 
Vo{n,M) vo + vo{n,M), 

where Jiq and Vq are the solution of the mean-field saddle-point equations 

_ h{ho) 

ho = 2{d - l)Pvl , 



(4.8) 
(4.9) 



(4.10) 



where c? = 5 in our case. If we insert the mean-field solution in the zero'th order 
effective action Eq. (14. 7p we obtain the zero'th order expression for the free energy 
per lattice site 

(4.11) 



• we have to fix the gauge; we take the Lorentz gauge which in the continuum amounts 
to adding the gauge fixing term l/(2^)9M^M9ArA^ to the Lagrangian. As men- 
tioned, the lattice version of this is Eq. (13. 6p . 

• the derivatives of the effective action are derivatives with respect to the components 
ha{n,M) and Va{n,M) . 

• the derivatives of the observables are derivatives with respect to the components 
4.2 The Faddeev-Popov determinant 



On the torus where the mean-field value of the links is universal, the ghost propagator is 
a simple application of Eq. (I3.12p : 



^ -I 5^^V {dl,d,A (4.12) 



M 



We then have (for SU{N)) the Faddeev-Popov term 



(det (wo^m^m)) 



Ar2-1 



In momentum space we will use the form 



^FP 



. p \ M 



Pll 



pM = 2sin(pAf/2) . 



(4.13) 



(4.14) 
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In the above, apart from a Fourier transformation we dropped a power of —i from App 
since it merely contribute an overall constant to the free energy. Instead, we choose to 
keep the factor of with the analogous consequences, as discussed in sect. 12.2.11 



4.3 The propagator 



Next we compute the matrix Eq. f l2.19p in Fourier space. The kernels which define the 
propagators in coordinate space are matrices 



K{n', M', a'] n" , M" , a") . 



(4.15) 



In Fourier space they are 

kip', M', a'; p", M", a") = ^Yl e'^'"'+<^'/'e-*^'"""-'^M"/2i^(n', M', a'; n", M", a") , 

n' ,n" 

(4.16) 

where the factor l/N" has been inserted to guarantee that the Fourier transform of the 
identity 6n'n" is the identity 5p'p" and tT{K} = tT{K}. 

The matrix K^'"'^^ is equal to the unit matrix. A straightforward calculation yields the 
result 

k{p',M',a';p",M",a") = Vp"<^c.'«"Ca/'m"(p', «') , (4-17) 

with the notation 



Cm'M"{p' ,Oi') — [A6m'M" + Bm'M"{^ — ^M'M")] , 

and (6i and 62 are defined in Eq. f l4.2ip below) 



A = - 

and 

Bm'M" 

where 



— (1 — Sa'o) + -;-Sa'0 

62 Ol 



cos [p'^) + (1 - (5a'o)^ sin^(pAf'/2) 



(4.18) 



(4.19) 



■ITT 2 



-4/3tJo 



c- I Pm'\ f p'm" 

Oa'o COS I ^f- 1 COS ' 



y/2(l-5.,o)sin(%^ sin 



(4.20) 



h. = - 



ho 



(4.21) 



and y = 2 — 1/^. For convenience we rescaled the gauge fixing parameter 1^. This 
propagator was presented in jl] in the axial gauge. 
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4.4 The Free Energy 



If we write Eq. fl2.16p in momentum space, the free energy Eq. f l2.28p can be evaluated as 



F(i) = F(0)+ 1 In 
2U 



JJdet 



Mhh) r^{vv) . , 



1 I i^-W i^-H) 



(4.22) 



The toron contribution p = is dropped according to our regularization scheme. 

The determinants in Eq. (14.221) are invariant under the flip of the sign of any of the 
components of the momentum p. We compute their contributions to the free energy by 
summing the logarithms of the absolute value of their eigenvalues and we keep track of 
their signs. The flip symmetry then implies that the overall sign of the determinant is 
decided by the momenta whose components are or vr. 

4.5 The Static Potential 



The observable we will use to compute the static potential is the Wilson loop extending 
ns points along the extra dimension and t points in time. We have defined the observable 
in eq. (I3.16p . The objective here is to obtain the correction as it appears in Eq. (12.3 
The calculation of Fig. [1] yields the first order correction to the tree level result 



0[V] = 2ivo) 



2{t+n5) 



(4.23) 



which can be brought in the form 



1 t 



2(- )2(t+ns)-2 J2 



P'l'P'2'P'i'P5'P'o=0 

[2 cos{p',n,) + 2]CoV(p', 0) + 3[2 cosip'.n,) - 2] 



Cooip',1] 



(4.24) 



Here we have used the fact that for a' ^ and p'q = the matrix K is block diagonal, in 
particular Cqq = I/Cqo- This can be seen by inspection of the matrix Ka'^o = 



/ -4 + ^'cAf - |sm"(^ 



-2(3v 



2(Po 
2 



ysi/2So/2 - 4 + X]' CM' - I sin 



yS2/2S0/2 yS2/2Sl/2 
ySz/2-%/2 yS3/2Sl/2 yS3/2S2/2 



ySo/2Sl/2 ySo/2S2/2 2/50/2^3/2 ySo/2S5/2 \ 

2(P[ 
5 -"^ I 2 



ySl/2S2/2 2/S1/2S3/2 2/S1/2S5/2 



|sin2(^) 2/S2/2S3/2 2/S2/2S5/2 



sm 



2{Ps- 



yS3/2S5/2 



jp'p" , 



\yS5/2S0/2 2/S5/2S1/2 yS5/2S2/2 yS5/2S3/2 - 4 + ^l' Cm' " 7 slu^ (^) / 



Y2' means sum over all values of M' except the one that corresponds to its column/row 
index (the index ordering is M = 0,1,2,3,5) and sm/2 = sin(p'^j,/2) (similarly, cm' = 
cos p'j^ and cm/2 = cos{p'^, /2)). 
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A crucial cancellation happens for a' 7^ 0. As can be seen directly from Eq. f l4.19p 
in this case, if = then A = and the matrix C has a zero mode. This "toron" is 
precisely cancelled by the factor [2cos{p'^n5) — 2], where the —2 comes from adding the 
self-energy and tadpole contributions. Another important observation is that the static 
potential which selects the value = is ^ independent. Also, it is easy to check that in 
the limit 

/3 — ^00, vo — ^ 1 (4.25) 

one recovers the proper form of the perturbative static potential. This is because the 
a' = part of the propagator is of order with respect to its a' ^ part (which itself 
is of order 1//3), for large (3. 

For the static potential we arrive at the final expression 



V{r) 



-21og(Uo) 

1 1 
2^1^ 



X 



-E 

4 ^ 



+3 



E 

^5^(2cos(p»-2) 



(2cos(p^r) + 2) 



•-^00 



ip',0) 



Coo(p',l) 



(4.26) 



In the above we have used translational invariance to average over Wilson loops of spatial 
size r along the M = 1, 2, 3, 5 directions. 



4.6 The Scalar and Gauge Boson masses 

In this section we compute the scalar and vector masses to the order illustrated in Fig. [2] 
and Fig. ^respectively. At this order there are no ghost contributions. The computation of 
the mass of the scalar observable on the torus starts from Eq. fl2.3ip . the contraction of the 
second derivative of the observable with the propagator. The observable is computed using 
Eq. fl3.2ip and Eq. f l3.20p Fourier transformed and the propagator given by Eq. fl4.17p . 

We define the very useful quantitjj^ 

A("^)(n5) = V ^ , f = r + 1/2 , (4.27) 
and its Fourier transform 

A(-3)(p)= ^ (4.28) 



^ The argument of iJo(w) refers to a more general situation where the background can depend on the 
position along the extra dimension, like on the orbifold [6j. Integer values for m label four dimensional 
links and half- integer values the links along the extra dimension. On the torus vo{r) = vos- 
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The final expression for tfie scalar mass observable in Fig. [2] turns out to be 

Ci\t) = ^{Pl^y Yl (Pot) E \^^'''\p,)fK'' {{p'o, 0,p',), 5, 0; (K, 0,p',), 5, o) , 

p'o P's 

(4.29) 

where Pq^^ is the Polyakov loop evaluated on the zero'th order background. It is interesting 
to note that the component of the propagator that contributes to Eq. fl4.29p in gauge space 
does not have a pole at the toron p' = 0. The toron contribution is in fact the dominant 
one. Also, because the only component of the propagator that gives a contribution to 
the mass is along the 0-direction and the latter is not gauge fixed, it is exactly gauge 
independent. 

For the W gauge boson mass we define the contraction 

K-\ip'o,f),5,a) = J2 A(^^)(K)A(^^)(-p'5)^''(/,5,a;j9',5,a) 

P's^P'^ 

(4.30) 

with the property 

K-\{p'oJ),5,a) =K-\{p'o,-if),5,a). (4.31) 
In addition, we define a second useful contraction as 

T~\tJ,a) = J2e''''^'K-\{p'oJ),5,a). (4.32) 

p'o 

A calculation similar to the one for the scalar, gives for Fig. |3]the result 

= ^(n^°^)'(^o(0))^X^5^sin2(p',) [W\tJ,l)f . (4.33) 

p' k 

The toron here is cancelled due to the sin^ (p'^) factor. Here gauge independence is not 
obvious because it is gauge fixed components of the propagator that contribute. A careful 
numerical investigation shows that the vector mass is actually gauge independent to an 
accuracy of 10^^^. 

In Appendix \M we give more details about the calculations leading to Eq. (14.29^ and 
Eq. Km . 

4.7 Anisotropy 

The only regime of parameters where dimensional reduction is guaranteed is when the 
physical size of the extra dimension becomes much smaller than the spatial dimensions, 
R « I. On the lattice, such a scenario can be realized when an anisotropy parameter 
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is introduced along the extra dimension. To this effect, following [15], we define the 
anisotropy factor 

7 = M (4.34) 

whose zero'th order value is 7 = a^/a^,, the ratio of the lattice spacings along four dimen- 
sional hyperplanes and along the extra dimension and 

2Na5 ^ 2Nal 

Pa = , /35 = ^ . 4.35 

Then, 

/35 = /37, Pa = -. (4.36) 

7 

A quantity that is worth keeping in mind is the four dimensional effective bare coupling 
which rescales to 

9l = ^ = ^-^ (4.37) 

27iR N^P ^ ' 

and it can in principle become non-perturbative for large enough anisotropy factor if the 

other parameters are kept fixed. We are also interested in the ratio 

iVs 27ri?a4 2t^R 

— = = —-7 (4.38) 

L I I 

which controls dimensional reduction via compactification, provided that it is kept fixed 
while 7 is increased. We distinguish three regimes of the anisotropy parameter when 
N,/L = 1: 

• 7 = 1. This defines the isotropic lattice which, in the limit of an infinite lattice, 
should represent the "non-compact phase " of the continuum gauge theory. 

• 7 > > 1 . This is a regime where the size of the extra dimension is small with respect 
to the spatial length and in the limit L — )■ 00 it represents the "compact phase" of 
the continuum gauge theory. 

• 7 < 1. This is a situation where in the L — )■ 00 limit the continuum gauge theory 
has a large extra dimension. If dimensional reduction is realized in this phase, it 
must be due to a localization effect. 

The 7 7^ 1 cases are interesting in case there is dimensional reduction. Then one ex- 
pects the four dimensional effective coupling to behave as a confining/ asymptotically free 
coupling. If this is indeed the case, such a property should be reflected by the static 
potential. 
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The parameters we choose to parametrize anisotropy are the couphng j3 and the 
anisotropy factor 7. In terms of these, we rewrite the Wilson action as 



^ 4d-p 5d-p 



(4.39) 



where the first term contains the effect of all plaquettes along the four dimensional slices 
of the five dimensional space and the second term contains the effect of plaquettes having 
two of their sides along the extra dimension. 

The partition function in the presence of anisotropy gets modified. All derivatives of 
plaquettes containing a link pointing in the extra dimension will be now different. There 
will be two mean values for the links, Vq and ^05 determined by the extrimization of 



S,s[VH]_ (d-l){d-2)_, 



— -P4- 



which yields the conditions 



VQ-l35{d-l)vQVQ^+{d-l)u{ho)+u{ho5) + {d-l)hoVo+ho5Vo5 

(4.40) 



— — 6/3 

Vq = -u{hoy , ho = + 2^7TJoVo5 , 



7 



^^05 = -u{ho5y : 



05 



8/37^0^05 , 



(4.41) 



for d — 5. 



The anisotropy factor cancels out in scalar products when written in terms of dimen- 
sionless momenta and coordinates 

pn = psHs + ^p^^n^ = (PsOs) (2^5/05) + J^(^m«4) (2^^/04) , (4.42) 

where in the second equality we explicitely put the lattice spacings. 

The modifications to the propagators are slightly more involved. The propagators in 
the anisotropic vacuum become 

/ 1 \ 



j^{hh) — — 5n'n"Srv'fy" = — ^ — — 



1-e 



hn 



hiho) Khiho) 



mho 



/o 



-dpipii da' a" = — .-r . 

"'05-'l("'05j 



1-e 



^'05 / -^2(^05) 



-^2(^05) V -^1(^05 



hihob) 



\ 



(4.43) 
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where e = 1 for a' = and e = for a' 7^ 0. For later reference, we note that the hmit of 
the factor in the second hne of the above expression for /iqs is —1/4. Also, 



tA'"'") c- c- / r,P —2\ 



( J2'cm' - 

ySl/2So/2 
I/S2/2S0/2 
I/S3/2S0/2 
Z/5-55/2S0/2 



le2 
^■'0/2 



yS0/2Sl/2 



where y 



Yl' CM' - 

yS2/2Sl/2 

I/S3/2S1/2 

1/5-S5/2S1/2 
1 



1,2 

ri/2 



2/ ■50/2^2/2 

ySl/2S2/2 

_ 1.2 
5 ''2/2 

I/S3/2S2/2 



2/S0/2S3/2 

2/ -51/253/2 

yS2/2S3/2 

is2 
5 ■'3/2 



2/5 



Z/5-S5/2-S2/2 Z/5S5/2-S3/2 7 
^^0 C 



2/5S0/2S5/2 
2/5S1/2S5/2 
1/5S2/2S5/2 
2/5S3/2S5/2 



lo2 
€■^5/2 



(4.44) 



and 



M^2, 

7' 



2Co/2Ci/2 2C0/2C2/2 2C0/2C3/2 2Co/2C5/27^ 

2lJ05 



2C1/2C5/27 

2C2/2C5/27'^ 



2Ci/2Co/2 Yl' ^M' 2C1/2C2/2 2C1/2C3/2 

2C2/2C0/2 2C2/2C1/2 Xl'cAf 2C2/2C3/2 ^/^ ^/^ , „Q 

2C3/2C0/2 2C3/2C1/2 2C3/2C2/2 E'cM' 2c3/2C5/27^ff 

V 2C5/2Co/27'|f 2c5/2Ci/27'|f 2C5/2C2/27' |f 2C5/2C3/27' |f 7^ E' ^M' / 

(4.45) 



The only point not explicitly shown in these formulas is that in the diagonal elements, 
C5 = 7^^ cos (p'5). 

In Eq. fl4.44p we have already implemented the anisotropic version of the gauge fixing 
term 

2 



^ A=l n L M 



'^{vA{n,ii) - VA{n - fi, fi)^ +-f(^VA{n,5) - VA^n - 5,5) 



(4.46) 



(up to a redefinition of ^). The resulting Faddeev-Popov determinant splits accordingly 
to two parts: 



A 



FP 



sm 



-) +7^05 sin 



2/«5P5- 



(4.47) 



As far as the observables are concerned, the modifications to the free energy and the 
mass formulae are simple. For the free energy Eq. fl4.22p all the information about the 
anisotropy is contained in the propagator and the Faddeev-Popov determinant. The cor- 
relators Eq. (14.291) and Eq. (I4.33p for the masses do not change. On an anisotropic lattice 
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there are two inequivalent Wilson loops. One along the four dimensional hyperplanes for 
which the static potential is given by 



V{r) 




+3 E 



^5^(2cos(p',r)-2) 



Coo(p',l) 



(4.48) 



The other Wilson loop is along the extra dimension and gives the potential 

V{r) = -21og(T;o) 
1 1 



I [2 cos(p;,r) + 2] Coo^(p', 0) + 3 [2 cos(p;,r) - 2] ^r^^^} • 



The matrix C is defined as in Eq. f l4.17p . 

For completeness and as a check, the perturbative limit /3 — )■ oo, — ^ 1, ^'os — )■ 1 of 
the static potential at distance r in direction M is easily computed to be 



cosip'^r) - 1 



«4 p/^^/=o ^kC^ ~ cos(p'fca4)) + 7^(1 - cos(p^a5)) 



(4.49) 



and for 7 — )■ oo one is left with a pure Coulomb potential at = (for clarity we have 
explicitly put in the lattice spacings). 



5 The phase diagram 

In the anisotropic background the phase diagram can be split in three regimes according 
to the solution to which the numerical iteration Eq. fl4.4ip converges. Where it converges 
to vq = = we define the confined phase. Where it gives a solution of the form 
^0)^05 7^ we define the deconfined phase and where it gives a solution of the form 
Vq 7^ 0, Uo5 = we define the "layered phase" pj. The latter is a regime where the four 
dimensional coupling is deconfined and the five dimensional coupling confined. This is 
an approach that gives a first approximation to the phase diagram. One then has to check 
whether the solution of Eq. fl4.4ip indeed corresponds to the global minimum of the free 
energy Eq. IK22\i . 

The mean-field method is blind to the confined phase of the five dimensional theory 
so we do not have to say anything about it. According to the first order free energy 
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P=0.905,Y=1/3.8 



10^ 




Figure 4: Left: The free energy at first order Eq. (j4.22p computed on isotropic T = L = = 10 
lattices as a function of /3 and the background value v. For /3 > 1.55 the minimum is at f > 0.80 
(at 0th order we found a non-trivial minimum for /3 > 2.12). For /3 < 1.55 the minimum jumps 
to values 0.42 <v< 0.48. Right: The free energy at first order in the layered phase for /3 = 0.905 
and 7 = 1/3.8 computed on a T = L = = 10 lattice. The minimum found at 0th order 
(uo = 0.9030, uo5 = 0) is stable. Where the free energy is complex, it is not plotted. 

P=2.136, 7=0.25 P=1,y=10 




Figure 5: Left: The free energy at first order Eq. (j4.22p in the d-compact phase for /3 = 2.136 
and 7 = 0.25 computed on a T = L = A^5 = 10 lattice as a function of the background values 
V and f5. The minimum found at 0th order (Uq = 0.9679 and vqs = 0.0288) is stable. Right: 
The free energy at first order in the compact phase for /? = 1 and 7 = 10 computed on a 
T = L = = 1{) lattice. The minimum found at 0th order {vq = 0.9179, wos = 0.9773) is not 
stable, it moves to a "phase" with = 1 (and v = 0.86). Where the free energy is complex, it 
is not plotted. 
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the layered phase is stable. The 0th order background (or saddle-point) solution vq = 
0.9030, ^05 = remains as a minimum when 1st order corrections are included, see right 
of Fig. 12 

For now, the deconfined phase will be the regime of our interest. The isotropic lattice 
is stable as the left of Fig. H] shows. Clearly this is a case not interesting from the point 
of view of dimensional reduction so we will not have to say much about it. The right plot 
of Fig. O shows the free energy for /3 = 1 and 7 = 10 (compact phase) computed on a 
T = L = A^5 = 10 lattice as a function of the background values v and ^5. The component 
of the 0th order background along the four dimensional hyperplanes vq = 0.9179 moves 
to 1; = 0.86 (precision 0.01). The background along the extra dimension however, which 
at zero'th order is V05 = 0.9773 is unstable since the minimum moves to V5 = 1.00 (with 
an uncertainty less than 0.01) when the 1st order correction is included. In fact, the same 
instability is observed everywhere inside the compact phase at small (3: the Wqs component 
of the background runs to 1, the value that is the solution to the background equations 
at /3 = 00. This instability could be an inherent property of the system or just a property 
of the first order approximation, we can not tell for sure until a higher order correction 
to the free energy is computed. 

The left plot of Fig. |5] shows the free energy for /3 = 2.136 and 7 = 0.25 computed on 
a T = L = A^5 = 10 lattice. As it will turn out this is in an especially interesting regime 
of the deconfined phase which we call the 'd-compact' phase. The 0th order background 
Vq = 0.9679, = 0.0288 is stable, including the first order correction we get a minimum 
at essentially the same values. Furthermore, this result is quantitatively unchanged on 
aT = L = A^5 = 6 lattice. In general, the d-compact phase, for 7 not too small, seems 
to have a deeper from the zero'th order minimum at f = and = 1. This global 
minimum however gradually disappears as 7 is lowered, as it is the case on the left of Fig. 
13 The local, first order minimum on the other hand is consistent with the zero'th order 
result and is separated from the global minimum by a regime where the free energy is 
complex. We conclude this part of the discussion by saying that all the caveats regarding 
the interpretation of the ^-dependent free energy mentioned before, remain. 

Next, using our local (Wilson loop) and global (Higgs and W masses) quantities we 
will characterize more precisely the various regimes of the phase diagram. We keep always 
A^5 = L. The first observable we use is the W gauge boson mass from Eq. ( 14.33p . In Fig. [6] 
we present a scaling study of the behavior of the mass in lattice units as a function of 
1/L for the point f3 = 2.136, 7 = 0.25. The errors represent the uncertainty in the 
plateau value of Eq. f l2.42p . A linear fit in 1/L gives an extrapolation a^my/ = 0.0006(3) 
for L — )■ 00 with x^/d.o.f. = 0.32. The slope is 12.51(1). Adding a quadratic term to 
the fit confirms the conclusion that the mass seems to be consistent with zero within 
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p=2.1 36, 7=0.25 




Figure 6: The gauge boson mass as a function of A linear extrapolation in the infinite L 
limit gives zero within errors. 

errors. We can not exclude of course a tiny (possibly exponentially suppressed) mass. 
We will come back to this issue when we analyze the static potential. A scan of the 
phase diagram reveals that the W mass is actually nearly independent of any parameter 
other than L. There is a very mild /3 and 7 dependence but it does not seem easy at the 
moment to attribute any quantitative physical significance to it. Therefore, we conclude 
that the deconfined phase, at least in regimes where the system is five dimensional (i.e. 
not dimensionally reduced), is in a Coulomb phase. Note that since the background is 
along the Hermitian component of the links it can not possibly play the role of a symmetry 
breaking vev. That is, in the absence of an explicit vev we can not probe the Higgs phase, 
unlike in a Monte Carlo simulation where choosing the appropriate input parameters they 
system can be in a Higgs phase or not. 

The second of our observables that we will use is the static potential derived from 
Wilson loops oriented along the four dimensional slices of the lattice which are orthogonal 
to the extra dimensioij^. The points obtained will be fitted to both four dimensional 
(co + c/r, where cq and c are the fitted parameters) and five dimensional {do + d/r^, where 
do and d are the fitted parameters) Coulomb forms. This will decide wether the theory is 
in a dimensionally reduced phase or not, as already mentioned. Let us fix /3 = 2.136. In 

^ Note that on the anisotropic lattice there are two types of inequivalent Wilson loops. One along the 
four dimensional slices and one along the extra dimension. We discuss here the former. 
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Figure 7: The static potential on the isotropic lattice. It is fitted by a five dimensional Coulomb 
form, as expected from the corresponding continuum theory. 



order to obtain a first, qualitative picture of dimensional reduction we perform global fits 
to four dimensional or five dimensional Coulomb potentials and assign to the potential 
points fictitious relative error^ and perform a comparative study of the of the fits. 
A much more careful study of the local curvature of the potential will be performed 
later. A scan in 7 of the global fits shows that for 0.35 < 7 < 4 there is no dimensional 
reduction. In particular on a lattice with L = = 100 for 7 = 1 one sees a perfect 
global five dimensional Coulomb fit, see Fig. [3 For 7 > 4, as expected, we start seeing a 
four dimensional Coulomb law which becomes better as 7 increases when only the points 
corresponding to small r are fitted. For large r we see deviations. We call this phase the 
compact phase. According to the free energy analysis the compact phase for small /3 is 
unstable to this order in the mean-field expansion, therefore we do not analyze further its 
properties. The somewhat surprising fact is that for 7 < 0.35 the potential turns again 
four dimensional Coulomb, until the layered phase is hit, at around 7 = 0.25. A similar 
picture can be obtained for other values of /3 in the deconfined phase. The narrow band 
for 7 < 1 which, at least according to our mean-field approach, apparently describes a 
phase where the extra dimension is large but the system is dimensionally reduced, we call 
the "d-compact phase". 

We now turn to the Higgs observable. By construction it is L-independent so we do not 
^ This is justified if we want to compare fits at different values of 7. 
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Isotropic case y=1 
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Figure 8: The scalar mass in lattice units on the isotropic lattice shows a ferromagnet type of 
behavior with exponent v = 1/4. 
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Figure 9: The scalar mass in lattice units a^^rriH in the d-compact phase shows a ferromagnet 
type of behavior with critical exponent v = 1/2. 



have to scale it with L. It turns out to be also independent and to have a strong /3 and 
7 dependence. At a generic point in the interior of the phase diagram the scalar is heavy 



32 



Y=0.25, p=1 .44333 
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Figure 10: Tuning of the bare coupling /3 as a function of L when the phase transition is 
approached at constant p = mw/rnH = 1.44333 and 7 = 0.25, = L. The sohd hne is a fit 
according to Eq. (j5.4p . 



in lattice units. It starts loosing weight as boundaries between phases are approached. 
In Fig. [8] we show its behavior as the isotropic phase transition is approached. The way 
the scalar mass approaches zero reminds us of the magnetization of a ferromagnet as a 
function of the temperature as the Curie temperature is approached where a second order 
phase transition occurs. Even though in our case the phase transition is definitely known 
to be of first order, the resemblance of the mass dependence is close enough so that we 
attempt to fit the data points to a critical law of the form 



arriH 



(5.1) 



which yields p = 0.2470 when we take (3^ = 1.6762017. When we take (3^ = 1.6762016 
we obtain u = 0.2546 instead. The average value from the two fits (0.2546 + 0.2470)/2 = 
0.251 ± 0.004 allows us to safely approximate the exponent to be 



1 

"=4 



(5.2) 



That this can not be a second order phase transition is also supported by the fact that the 
mass does never really go to zero. It reaches a small value and it stops there. The same 
behavior can be verified along the whole line that separates the confined and deconfined 
phases. 
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A qualitative difference appears on the boundary between the d-compact and layered 
phases. To begin, the scalar mass in lattice units goes all the way down to zero as the 
layered phase is approached. Moreover, in Fig. [9] we measure a critical exponent 

^ = \ (5.3) 

when we set /3c = 2.1349, indicating that the d-compact phase describes different physics 
from the rest of the deconfined phase. We note that this is the critical exponent of the 
four dimensional Ising model and it is believed that if the SU{N) theory should have a 
second order phase transition in d+1 dimensions, its order parameter is the Zj^r values 
of the Polyakov loop and its critical exponent should he u = 1/N \T7]. Furthermore, if 
there is a tricritical point then above it the phase transition is expected to turn into a 
first order one with "exponent" 1/4. We have also checked that the free energy, as the 
phase transition is approached looses its local maximum and it becomes flatter around its 
global minimum, both facts consistent with a second order phase transition [5]. In our 
case the second derivative of the free energy at the minimum approaches zero as the phase 
transition is reached. Whether this is the real physical picture or it is just an artificial 
effect that the mean-field method produces, is not clear. Until a Monte Carlo simulation 
study is performed we can not be sure. Finally, we have examined the passing from the 
layered phase into the confined phase. As the scalar masj^ gets stuck at a large value, 
this can not be a second order phase transition. In this case the indication is that it is a 
first order phase transition, like in the U{1) theory [T8] . 

For now we assume that what we observe is a second order phase transition and we 
study the scaling behavior as we approach it. First, we tune the values of /3 as the phase 
transition is approached keeping p = mw/^H fixed (we keep = L and 7 fixed as 
usual) as a function of the lattice size, see Fig. [TOl § The behavior fits well to the form 

/3 = Co + -^ = co + ci^, (5.4) 

which implies that the bare couphng squared gl = A/ [3 approaches its critical value as a 
quadratic function of the lattice spacing. 

Second, for the values of I3{L) at constant p = 1.44333 we calculate the static potential 
V{r) in the four-dimensional slices orthogonal to the extra dimension and fit it to the four 
dimensional Coulomb form 

V(r) = h const . (5.5) 

r 



^ The Polyakov loop Eq. p.lSp is zero in the layered phase. In order to measure the scalar mass there, 
we rotate the Polyakov loop in the k — 3 direction. 

^It is important to notice that the ratio p is independent on since both mw and rriH turn out to 
be independent on N5. This means that as the lattice spacing goes to zero as the phase transition is 
approached one can keep R fixed and therefore, when L — N5, also / fixed. 
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Figure 11: Scaling of the charge defined in Eq. ()5.5p as the phase transition is approached 
along the values /3(L) of Fig. [10] at constant p = 1.44333. 



The fits become better as L increases. The fitted parameter can be interpreted as a 
charge. Since the static potential is renormalized up to an additive divergent constant [T9] . 

is a physical quantity and Fig. [11] shows its scaling behavior. The points L = 16 ... 48 
closest to the phase transition can be very well fitted by 

5 0.30 , , 

q^ = — + 0.0342 5.6 

with x^/d.o.f = 3.7/4. The fit is shown in Fig. [11] Adding a 1/L term to the fits give 
g2 = _o.37/L2 + 0.0054/L + 0.0341 with xVd-o.f = 3.0/30 so we conclude that the 
continuum limit of is approached up to quadratic cut-off effects as expected [20] . 



5.1 The static potential 

We now turn to a closer examination of the static potential for 7 7^ 1. We concentrate on 
the d-compact phase and we choose to approach the 2nd order phase transition keeping 

7 = 0.25 

fixed. We choose the values /3 = 2.13495, 2.136 and 2.3, which correspond to Higgs masses 
of a4mH = 0.020, 0.094 and 1.12 respectively. In the d-compact phase, the relations (for 

^'^ We assign artificial relative errors of 0.001 to the potential in order to compare the fits. 
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Potential along the extra dimension, p=2.136, y=0.25 
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Figure 12: The static potential derived from the Wilson loop along the extra dimension is 
independent on and has a strong L dependence. As L increases, its slope decreases. 



N, = L) 

a5 = a4/'j and R/l = l/{2iT'y) (5.7) 

indicate that the extra dimension is large compared to the size of the four-dimensional 
space. We take the continuum limit at fixed anisotropy, as indicated by the arrow in 
Fig. [21] It is interesting that only if 7 < 1 a continuum limit exists, that means that 
the anisotropy survives in the limit and has physical consequences. We measure the 
static potentials along four- dimensional hyperplanes orthogonal to the extra dimension 
and along the extra dimension. 

We start by looking at the potential along the extra dimension. In Fig. [12] we plot the 
static potential at /3 = 2.136 for various values of A^s and L equal to 32 or 64. Due to 
the periodicity, we measure the static potential up to distances L/2 and we can therefore 
compare distances up to 16 in lattice units. One can see that it is essentially independent^ 
on A^5 and that its slope strongly decreases as L increases. These two facts give rise to 
the suspicion that it may not represent a physical interaction. To check this, we plot in 
Fig. [13 the physical dimensionless force 

f^Fif ) = [V{r) - V{r - ag)] /ag, f = r - ^a^ + O(a^) (5.8) 
derived from the potential, keeping = L and increasing L = 24, 32, 48, 64, 100 (left) and 
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Barring small effects due to the periodicity of the lattice. 
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Force along the extra dimension 
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Figure 13: The dimensionless force between two static charges separated in the extra dimension. 
Apparently it is a pure finite size effect that approaches zero in the infinite volume limit. 
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Figure 14: The static potential derived from four dimensional Wilson loops is and L inde- 
pendent. 



L = 100, 125, 150, 200 (right). The plots show that the force goes to zero as the L — )• oo is 
approached suggesting that the four dimensional hyperplanes decouple in this limit and 
the system reduces to an array of non-interacting "branes" where gauge interactions are 
localized. 

We continue by looking at the potentials along four- dimensional hyperplanes. First, in 
Fig. [Hlwe check and L dependence. As the figure shows, these potentials are essentially 
independent of and L. Thus, it is likely that they reflect physical interactions. Recall 
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Potential along the 4d hyperplanes, y=0.25, L=200 
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Figure 15: The static potential along the four dimensional planes computed on L = 200 lattices 
approaching the phase transition in the d-compact phase and keeping 7 = 0.25 fixed. 



that the global fits of the potential along the four dimensional hyperplanes and the analysis 
of the extra dimensional force suggested that in the large L limit the potential along the 
four dimensional planes should describe four dimensional physics to a good approximation. 
In order to show this we analyze the potentials on L = 200 lattices for /3 = 2.13495, 
/3 = 2.136 and /3 = 2.3, shown in Fig. [151 We perform local fits to the form 



c d 

V(r) = Co + ar -\ h^ + / ln(r) . 



The inclusion of a logarithmic correction is motivated for example in 
dimensionless potential a4V"(x), x = r/a^, we derive the auxiliary quantities 



Zi{x) = x^[V{x + l) -V{x-l)]/2 
Z2{x) 
Y,{x) 
Yi{x) 



Zi{x + 1) + Zi{x-1) -2Zi{x) , 
x'[V{x + 1) + V{x - 1) - 2V{x)] , i 
x''[Y2ix + l)-Y2ix-l)]/2, 



2,3,4, 



in terms of which we can estimate the coefficients in Eq. (15.91) to be 



a{x + 0.5) = a{x + 0.5)al 
c{x + 0.5) 
d{x) = d{x)/ai 
l{x) = l{x)a4 



[Z2(x + l)-Z2(x)]/6, 
-[Y^{x + l)-Y^{x)]/2, 
x^[Ys{x + 1) + Y^ix - 1) - 2y:3(x)]/12 , 
-[Y^{x + 1) + Yi{x - 1) - 2r4(a;)]/2 . 



(5.9) 
From the 

(5.10) 
(5.11) 
(5.12) 
(5.13) 

(5.14) 
(5.15) 
(5.16) 
(5.17) 
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Potential in the 4d hyperplanes, y=0.25, L=200 Potentiai in the 4cl hyperplanes, y=0.25, L=200 
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Figure 16: Analysis of the static potential along four-dimensional hyperplanes in the d-compact 
phase on lattices of size L = 200 at fixed anisotropy 7 = 0.25. Left: The coefficient c(r) of a 
four dimensional Coulomb term. Right: The coefficient d{r)/ai of a five dimensional Coulomb 
term. These coefficients are determined from local fits to the form in Eq. (j5.9p . 
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Figure 17: Analysis of the static potential along four-dimensional hyperplanes in the d-compact 
phase on lattices of size L = 200 at fixed anisotropy 7 = 0.25. Left: The string tension a{r)a\. 
Right: The coefficient l{r)ai of a logarithm term. These coefficients are determined from local 
fits to the form in Eq. (|5.9p . 



In Fig. [in] we plot the coefficients c and d of the four dimensional and five dimensional 
Coulomb term respectively. Away from the phase transition, at /3 = 2.3, we observe 
almost plateaus for the coefficients d G [0.15,0.18] and c G [0.004,0.007] in the range of 
distances x = 40 ... 60. As can be seen in Fig. [171 the coefficients a and I are basically 
zero in this same range of distances. Therefore, the potential is a balance between a 
four- dimensional and a five-dimensional Coulomb term and we interpret this as the onset 
of dimensional reduction from five to four dimensions. 

In fact, as we go closer to the phase transition at /3 = 2.136 and /3 = 2.13495 the 



39 
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Figure 18: Coulomb fits of tlie static potential along four-dimensional hyperplanes close to the 
phase transition in the d-compact phase. The range r/a^ = 10 ... 50 is fitted to a four of five 
dimensional Coulomb form. 

potential becomes four dimensional. In the range of distances x = 15 ... 35 we see in 
Fig. [TBI a plateau for the four dimensional Coulomb coefficient c G [—0.031, —0.029], where 
the coefficient of the five-dimensional Coulomb term is much smaller d G [—0.01, —0.005]. 
In the same range, see Fig. [TTl the coefficient of the logarithmic term I is basically zero and 
the string tension has a plateau at a very small but negative value a G [—3 x 10~^, —2 x 



In Fig. [18] we fit the potential at /3 = 2.136 in the range x = 10 ... 50 globally to purely 
four or five dimensional Coulomb forms (i.e. we set o" = (i = / = in Eq. fl5.9p ). Clearly, 
the four dimensional Coulomb is an excellent fit while the five dimensional Coulomb law 
is exlcuded. Also, these plots exclude a possible large contribution from the presence of a 
light excited vector state which would contribute an extra Yukawa term in the short (and 
therefore also in the long) distance part of the potential. As already briefiy mentioned 
in the beginning of this section, this was our first indication for dimensional reduction in 
the d-compact phase. Together with the vanishing of the extra dimensional force, it can 
be taken as a strong evidence for dimensional reduction via localization. 



10~% 



40 



5.2 Discussion 



We have now to compare these results with the behavior that we expect based on the 
underlying physical picture. If conventional Kaluza-Klein dimensional reduction works 
and the low-energy theory contains a massless gauge boson and a Higgs particle (the 
assumption here is that higher excitations are separated by a large mass gap), then this 
low-energy theory corresponds to the four- dimensional Georgi-Glashow {SU{2) gauge the- 
ory with adjoint Higgs) model. In [2^ the phase diagram of this model is discussed, it 
has a confined phase and a Higgs phase where the gauge symmetry is broken to U{1) [22] . 
Note that in this model string breaking cannot happen when the static charges are in 
the fundamental representation and the matter in the adjoint. In order to attempt an 
interpretation of the long distance part of the static potential that is consistent with such 
a physical picture we will assume that the data in this regime is not unphysical. Since 
we have no guarantee that this is the case, the following discussion must be read with 
great caution. Also, in the following we will take seriously the short distance analysis of 
this section which lead to the conclusion that the system in the d-compact phase is truly 
dimensionally reduced. 

On one hand, the fact that we do not see a non-vanishing mass for the gauge bosons 
(no Higgs mechanism) implies that our four dimensional effective theory should land in 
a confined phase. On the other hand one would expect a mass gap to appear in the con- 
fined phase which means that the mass of a gauge invariant vector operator like Eq. ( 13.24p , 
representing a vector glueball state, should be non-zero. To reconcile this apparent contra- 
diction we recall the argument of the first reference of [3] where the dimensional reduction 
of a five dimensional gauge theory 1^ in a Coulomb phase is discussed. There, in the di- 
mensionally reduced state, the non-perturbatively generated mass gap behaves like 



{k is a constant) when ^4 — )■ and a — )■ 0. We have denoted the asymptotically free 
coupling by to distinguish it from our effective coupling defined in Eq. fl4.37p and 
which approaches a non-zero critical value as a — )■ 0. The mass gap in lattice units 
is exponentially suppressed near the continuum limit. Thus, in order to interpret our 
results in the context of dimensional reduction we must assume that also in our case such 
a mechanism is at work. This is an assumption since our data is not able to distinguish 
a zero from an exponentially small mass. 

Based on these arguments we expect to see in the large distance behavior of the 
potential a string tension contribution, which rises linearly with the distance. Instead 

In the context of D-theory the scalar field is decoupled, this can be achieved by choosing orbifold 
boundary conditions. 
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Figure 19: Analysis of the static potential along four-dimensional hyperpianes in the d-compact 
phase on lattices of size L = 200 at fixed inverse correlation length aumn = 0.1. Left: The 
coefficient c(r) of a four dimensional Coulomb term. Right: The coefficient d{r)/ai of a five 
dimensional Coulomb term. These coefficients are determined from local fits to the form in 
Eq. (fOD . 



the local fits based on the assumed form of the potential give a very small but negative 
string tension, see Fig. [171 This could signal a breakdown of the meanfield at large 
distances, absence of dimensional reduction (these go against our working assumptions 
though) or that there are important terms in the confining string potential which we have 
neglected. In order to circumvent this ambiguity we have assumed the presence of the 
universal Liischer term — vr/ (8r) in c? = 5 dimensions |23], subtracted it from the data and 
performed global fits of only the long distance part of the rest, varying L. In this case we 
obtained a positive string tension. The same analysis applied to the long distance part 
of the isotropic potential shown in Fig. [7] gave a similar result though, with almost the 
same value for the string tension. This means that the ambiguity persists. 

In order to check whether it is possible to find a unambiguous signal for a positive 
string tension we computed the potential along the line of phase transitions separating 
the d-compact from the layered phase. We choose three points in the (/3,7) parameter 
space, represented by blobs in Fig.EU (1.358, 0.5), (2.136, 0.25) and (2.998, 0.172). These 
points are at the "same" distance from the phase transition if we measure it by the inverse 
correlation length a^rriH (the Higgs mass in lattice units), which for the chosen points is 
0.101, 0.094 and 0.103 respectively. We present the results for the coefficients c, d in 
Fig. [19] and for the coefficients a, I in Fig. [201 While not much difference is seen between 
the points (2.136, 0.25) and (2.998, 0.172), for the point (1.358, 0.5) closest to the tncntical 
point of the phase diagram (the point where the compact, layered and d-compact phases 
meet) we get a positive string tension in the range of distances r/a^ = 27... 43 with 
average value 1.2 x 10~^. The value of the string tension increases as L increases. We also 
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Figure 20: Analysis of the static potential along four-dimensional hyperplanes in the d-compact 
phase on lattices of size L = 200 at fixed inverse correlation length a^mn = 0.1. Left: The 
string tension a{r)a1. Right: The coefficient /(r)a4 of a logarithm term. These coefficients are 
determined from local fits to the form in Eq. (15. 9p . 



note that the coefficient c of the 1/r term increases by almost an order of magnitude as 
we get closer to the tricritical point. These observations deserve to be studied further. 

We conclude this section by summarizing our results for the phase diagram of the 
anisotropic five dimensional SU (2) gauge theory in Fig. [2T] and by making a couple of 
remarks regarding the validity of our mean-field results. 

Firstly, the results we presented for our physical quantities (static potential, masses) 
are independent of the gauge fixing parameter ^. This is an immediate consequence of 
the fact that these quantities are gauge invariant and provides a check of our programs. 
The free energy depends^ on ^ but not the conclusions we derive from it. We checked 
that the instability of the compact phase and the stability of the d-compact phase are 
independent on the choice of ^. The precise location of the minimum of F^^'^ depends on 
^: for example at /3 = 2.136 and 7 = 0.25 on a T = L = Ai's = 10 lattice, the minimum is 
at V = Vq and = 0.48 = 1) or v = vq and = 0.42 = 10). Thus we would not 
base a serious criticism of our results on that. 

Secondly, it is difficult to predict the effects of higher order corrections and/or the 
general domain of validity of the mean-field approximation. The mean-field is known to 
produce fake physics in certain cases, so it may generate consistent but unphysical phase 
transitions. This criticism concerns mainly the line of second order phase transitions 
separating the d-compact from the layered phase, where we have carried out most of our 
analysis. This is indeed possible but our working assumption has been that whenever 



The gauge dependence of the free energy is inherent in the meanfield approximation. Already at 0th 
order its value depends whether we fix or not the gauge. 
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Figure 21: The phase diagram for lattice SU (2) gauge theory on an anisotropic torus, according 
to the mean-field method. In bold, the line of second order phase transition of critical exponent 
1/2 is plotted. The bold dashed line that separates the layered phase from the confined phase is 
a line of first order phase transition. The normal line indicates also first order phase transitions. 
The dashed lines indicate cross-over regions. The gray shaded area is unstable according to the 
first order free energy. The arrow shows the approach to the continuum limit for the discussion 
of the static potential at fixed anisotropy in Sect. 15. li The blobs represents three points at fixed 
inverse correlation length a/^mH discussed in Sect. 15.21 

the mean-field is non-trivial, it describes a physical property. The existence of such an 
ultraviolet fixed point is suggested by the epsilon expansion [27] but has been elusive so 
far in Monte Carlo simulations. This also deserves further study. 



6 Conclusions 

Using the mean- field approach we explored the phase diagram of a five dimensional SU{2) 
gauge theory with periodic boundary conditions. On the anisotropic lattice, we found two 
regimes where dimensional reduction seems to be at work. 

One, the compact phase at large values of the anisotropy parameter is found to be 
unstable for moderately small values of (3 down to the phase transition into the confined 
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phase when the free energy is computed at first order. In the hmit /3 — t- oo the mean- 
field reduces to standard perturbation theory where the compact phase does not have an 
instabihty but the theory becomes trivial. 

The other regime, the d-compact phase, is at small values of the anisotropy parameter 
as approaching the phase where the coupling along the fifth dimension becomes confined 
- the layered phase. The latter phase is metastable to first order in the free energy. 
The transition between the d-compact and the layered phase is found to be a second 
order phase transition. This allows us to take the continuum limit of physical quantities. 
We explicitly computed the continuum limit of a physical charge. In a dimensionally 
reduced state we are entitled to view the spectrum in a basis of four dimensional quantum 
numbers. In this basis, ours is an SU{2) gauge-Higgs system with the latter in the adjoint 
representation. Approaching the critical line, the Higgs mass remains finite while we are 
removing the cut-off. This is the first example of a four dimensional theory with this 
property as far as we know and is an independent check of results from Monte Carlo 
simulations [25] and effective field theory approaches [26]. The gauge boson mass vanishes 
within errors in the limit of infinite volume at any fixed lattice spacing. We cannot exclude 
an exponentially small mass expected if the system reduces to an effective four dimensional 
theory in a confined phase. 

We investigated dimensional reduction in the d-compact phase by computing the static 
potential. The static potential along the extra dimension seems to vanish in the infinite 
volume limit, a fact which would support a scenario of dimensional reduction via local- 
ization in the d-compact phase. Moreover we found that as we get closer to the phase 
transition, the potential along the four dimensional hyperplanes turns into a four dimen- 
sional Coulomb form at short to intermediate distances. The unambiguous observation 
of a positive string tension contribution to the mean-field potential along the four dimen- 
sional planes at large distances remains elusive. But as the tricitical point is approached, 
local fits indicate a small but clear plateau were the string tension is positive. 

It remains to see what a fully non-perturbative Monte Carlo computation will tell, in 
particular concerning the existence of a second order phase transition at small values of 
the anisotropy parameter. Our next analytical step following this work is to change the 
boundary conditions along the fifth dimension from periodic to orbifold |16] . 
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A Meanfield corrections to Scalar and Gauge Boson masses 

In this Appendix we describe in detail the computations of the diagram in Fig. |2] and 
Fig. [3] leading to the formulae given by Eq. f l4.29p and Eq. fl4.33p respectively. 



A.l The Scalar mass to first order 



The first order correction to the Higgs mass is given by Eq. (I2.3ip with the operator O 
defined by Eq. fl3.20p . The second derivative to be computed is 

5tT{p(^\to,m')}5tT{p(^Kto + t,m")} _ (^,) (^,) (o).2 

where Pq^^ is the background value of the Polyakov loop and the symbol A^^'^^n^) is 
defined in Eq. fl4.27p . After averaging over starting time to and space positions m', m", 
the Fourier transform is 

ni,ra2 to rh' ,m" 

^^(no)ify+tS(no)2,to^ni,rn''Sn2,rh' + S{no)iM^{no)2,to+t^ni,m'Sn2,rh'' 

to rn' ,m" 



to 



(A.2) 



We have also added the term that corresponds to interchanging the index 1^2. Putting 
everything together and contracting with the propagator, we have the final expression for 
the first order scalar mass observable Eq. f l4.29p . 



A.2 The Gauge Boson mass to second order 

Next we turn to the W bosons. As discussed, the contribution to the gauge boson's 
masses comes only at second order and is given in Eq. f l2.43p with the operator O defined 
by Eq. fl3.24p . What we will need is the derivative of the displaced Polyakov loop wjf^'"^ 
defined in Eq. fl3.22p . which is essentially determined by the single derivative of 
Eq. fl3.19p . All terms with derivatives acting on the links U and any terms with more 
than one derivative acting on vanish identically because of = when evaluated 
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on the background. The single derivative is simply 

2S"Po^°^'^„„i5„^,^5M5A(^^)K) . (A.3) 



6va{n,M) 

On the torus, we have to take S*^ = and S"^ = icr^0 The non- vanishing contributions 
to the fourth derivative in Eq. fl2.43p are of the form 

,5 iyi°^'^(to, m') S iy/°^-^(to + t, m") 

6va^ (ni, Mi)6va2{n2, M2) Svasim, M3)6Va^{n4^, Mi) 

A(^^)(K)i)A(^^)(K)2)A(^^)(K)3)A(^^nK)4) 

+ (1 2) + (3 4) + (1 2, 3 O 4) j , (A.4) 

where (1 ^ 2) means that (ai, ni) in the second last line of Eq. (IA.4p has to be exchanged 
with (02, ''^2) and similarly for the other exchanges of indices. According to Eq. (12.431) we 
have to compute the schematic expression 

1 V- 1 V- 1 V- „ 



(ni,Mi,ai)---(n4,M4,a4) to m' ,m" 

K~\l; 2)K'\3; 4) + K-\l; 3)K-\2; 4) + 7^-^(1; 4)i^-^(2; 3)1 (A.5) 

in momentum space. The factor 6 is a symmetry factor that accounts for the different 
choices of pairs of indices out of vi, f 2, fs and ^4 when taking the double derivatives. The 
first term in the bracket corresponds to time independent (self energy) contributions so it 
can be dropped, since it cancels out from the connected correlator. We will rewrite each 
of the two remaining terms as 

''^Oi{m;n)K^^{n;m) x 02{r] s)K~^{s; r) 

m.n r.s 

= tiiOiK-^) X tr(02i^"^) = tT{diK-^) X tr(02i^~^) • (A.6) 

For easier reference, we call the two relevant contractions 1 and 2 (the second and third 
term in the bracket of Eq. (lA.Sp respectively). Also, we label the four terms in the last 
part of eq. (lA.4p with labels from I to IV. 

In order to keep our formulae as general as possible we will not assume momentum 
conservation in the fifth direction (like on the orbifold) and we will consider a propagator 
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By introducing the matrices E", Eq. (|A.3p is applicable to the orbifold case as well. 
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that can have off-diagonal terms in gauge space (hke it can happen when the mean-field 
background is not simply proportional to the unit matrix). 

A term of the form 

0(m;n) = 0{mQ,nQ) O {171^,71^)0 {rn,n) (A. 7) 

will transform in momentum space into 

0{p';p") = d{p'„p'^^ e^"»"^°e-^^onoo(^^^^^)i_^e^p'■rn^-.p-^nQ^^^^ ^ 

mo, no rn,n 



(and analogously using q'; q" for the operator in the other factor of Eq. (1A.6P ) where the 
part along the extra dimension is the universal term 

Oip'^p'^) = ^AW^^)(p'5)AW^^)(V5') . (A.9) 
The temporal Fourier transform gives for the first contraction 

0(3; 1): le*Po(*»+*)e-^J'o*% 0(4; 2): Le<('o+t)^-iq'^to ^ (^^^0) 
and for the second contraction 

0(4; 1) : —e^Po(^o+t)^-iPoto^ 0(3; 2): iLe*'?o(*o+*)e^*'?o*o _ 

Both products result into 
to 

where we made conservation of temporal momentum, enforced upon taking the trace with 
the propagator, explicit. Next we turn to the spatial Fourier transforms. Each of the two 
propagator contractions contains four terms, which we list here inserting conservation of 
spatial momentum. 

Contraction 11: 

^ \ ^ ^ X X iff (m" +i—(m' +k)) i<f (m" —m') r e e „i(p',—p'i) 



Contraction III: 



Contraction II H: 



J2 ^(5p.,p..(5g.,,-».e*P"(™''-("'+^))e*9"«'^''+')-"') = ^6f,j?>Sg^,^,5j?,^,^e-'^'''^^P'^^ 
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Contraction IIV: 



J6 lE^f'f'^' 
rh'm" 



ip'{rn''-rn') iq'{{m''+l)-{rn'+k)) 



Contraction 21: 



L6 L6 



Contraction 2/J: 



1_ V — r^^ ^,r^^ ip'(m"-m')p»9'(("i"+0~(m'+fc)) 



L6 



Contraction 2///: 



L6 L6 



Contraction 2/1^: 



ip* ((m" +0 - (m' +fe)) gig* (m" -m' ) 



^'fe-' 



^g j^e"'p',p""q',q"^ ^ ^gWp',p"W(f ,9"^,-?"^ 



L6 



We can simplify the expressions if we define 

Z-'((p;,p-'),5,ai,a2) = 5^ A(^^)(p'5)A(^^)(-p^0^-'(P^5,ai;p',5,a2). (A.12) 



The components of the propagator diagonal in the Euclidean index are even under re- 
flecting the spatial momenta, which means 



(A.13) 



If we define in addition 



K {tj,ai,a2) = ^e'P'»'K ((p'o,p'), 5, ai, ^2) 



(A. 14) 



and use the fact that S"''" = — S° we can write out the whole expression Eq. ( 1A.4I) as 
6- 16-4 1 
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(Po^°^)^(l;o(0))^^ trja^S^^S^^jtrla^S^^S""! 



-1 



COS {pi-p^)K {t,p,ai,a3)K {t,p, 02,04) 



+ cos {p[+p'^)K {t,p',a2,a3)K 01,04) 



(A.15) 
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Finally we simplify this expression in some special cases. If the propagator is diagonal in 
the gauge index, Eq. (]A.15|) becomes 
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(n^°^)'(^o(0)r$^5^tr{a^S-S-}- 



'p' OL\^OL2 



ill S"^} ^cos (pO cos {p',) + tr^ S"^] ^sin {p[) sin {p[ 



K {t,p',ai)K {t,p',a2) 



(A.16) 



On the torus, S'^ = and T,^ 



imply that the term with the cosines in Eq. (IA.16P 



is absent and so the toron does not contribute to the correlator. Using ti{a^a^^a^'^} = 
2^g^AiA2 and summing over the gauge index A and the Euclidean indices with S^i gives 
the result in Eq. fOSj) . 
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